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ABSTRACT 


The Pacific Plate Motion Experiment (PPME.) will use intii ometry tech- 
niques to make precise measurements of geophysical v e- s . At the level of 
accuracy required, the path length error due to the neuUtu iov»er atmosphere 
makes a significant contribution to the error budget of the interferometer 
measurements. In order to better define the needs of the PPME, this review 
and studv was prepared. 

A review of previous work shows that, while very little has been published on 
the small-scale distribution of the atmospheric water vapor, the large-scale 
distribution of watex* vapor is known to vary significantly over time scales as 
short as 20 minutes. This variation in water vapor content is the main source 
of the changes fn path length error at a particular site. 

It is shown that the current "state of the art" in passive microwave radiometry 
allows a very precise measurement of the brightness temperature of the sky. 

It is also noted that the technological requirements of radiometers are very- 
different from the requirements of radio astronomy. The technology has been 
used in the construction of radiometers which are sufficient for use in the path 
length correction problem. 


A simulation study shows that, when combined with sui’face meteorology data, 
passive microwave radiometer data will allow' a determination of the path 
length correction to better than 2 cm at the zenith. By a careful choice of 
frequencies, a dual frequency system will allow a measurement of the path 
length correction to better than 4 cm at zenith angles as great as GU°. Be- 
cause of the wide range of weather conditions to be expected for the PPME 
sites (which Include Alaska, Hawaii and Massachusetts), it will probably be 
necessary to use a separate correction algorithm for each site. 
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I. THE PATH LENGTH PROBLEM 


The* atmospheric path length differs from the geometric |>ath length because of 
the influence of the index of refraction of the atmosphere. If 7 is the vector 
from the observer to the object of interest, the difference between the geo- 
metric and apparent path length t7T> is given by: 

b 

K« f m l)dr 

a 

where the observer is at point a, the object is at point b and n is the index of 
refraction. 

Any precise measurement of a length within the atmosphere requires a correc- 
tion for the path length error due to the atmosphere. Obviously, the impor- 
tance of the correction increases as the geometric path length increases until 
the ent r«* atmosphere is included For example, the path length correction 
for the zenith direction through the entire atmosphere is about 2.2 meters for 
frequencies around 3 GHz. 

In this paper, we will deal with the problem of path length corrections in the 
1.4 GHz to 10 GHz frequency range. As contrasted with the optical problem or 
the millimeter wavelength problem, the following conditions hold: 

1. DisjKjrsive (i.e. , frequency dependent) effects are negligible. Calcu- 
lations show that at most 10 '’of the |>ath length correction at 10 GH' 
is due to dispersive effects. 

2. Clouds are nearly transparent. Except for the highest frequencies, 
even heavy rain clouds are not opaque. 

.3. The dominant influence on the variation of the path length is the dis- 
tribution of atmospheric water vapor and its changes. As we will see, 
the index of refraction of the dry atmosphere is relatively constant, 
while the index of refraction due to water vapor is highly variable. 

In a sense, the dry atmosphere sets the average path length error 
while the water vapor is responsible for the variations. 

IU*an and Dutton (1966) have shown that the radio refractive index can be ap- 
proximated by: 
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where e is the partial pressure of water vapor (mb), l Is the lemiieruture 
( K) and P is the atmospheric pressure (mb). A and B are constants (A 77.fi 
°K 

x It) ’ B 4M0 ' K). Substituting into (1) ami using the ideal hub law 

mb 

yields, for the zenith direction: 

m 

Klein) = 0.277 ho) + I dh (3) 

) 1 2 ( h ) 

o 

where h is the altitude above the surface (cm) and C is a constant. 

Clearly, at a minimum, the surface pressure and the altitude distributions oi 
water vapor and temperature are required to determine It. 

Of tlie approximately 2.2 m of typical path length, equation (3) shows that 
about I. 5 meter is due to the d"’* atmosphere. Since seasonal variations of 
the partial pressure of water vapor can be as much as 100'* there can be vari- 
ations of nearly half a meter in tl e water vapor contribution to the correction. 

Numerous studies have obtained approximations to the path length correction 
based on atmospheric statistics and surtace meteorological measurements 
(i.e. , llopfield, 1971; Marini and Murray, 1073). I he conclusions have usu- 
ally been that the correction schemes are accurate to about 2 cm for the 
zenith. However, a recent study (Goad, 19741 shows that, when the various 
formulations are used tor a site which was not included in the original statis- 
tical ensemble, the predictions are accurate to 5 cm at best for the zenith and 
become progressively poorer as the zenith angle increases. Since the contri- 
bution of the atmospheric path length error to the accuracy of the Pl’MF base- 
line measurements cun be considerable, it is necessary to measure the ;«ith 
length correction to 2 cm or better accuracy at the zenith. 

Because of the difficulty of implementing active sounding techniques near a 
high sensitivity receiver, rada r -type systems cannot be used to measure the 
path length correction. Because clouds are opaque to visible and infrared ra- 
diation, optical sounders will not measure the entire geometric path except for 
clear skv conditions. We are therefore restricted to passive microwave tech- 
niques to evaluate the path length correction. 

The microwave spectrum of the earth's atmosphere contains two spectral line 
complexes oi importance for the path length problem. A measure of the dis- 
tribution <>f water vapor (and thus e(h) ) can be obtained from observations 
around the 22.235 C.llz water vapor line. Similarly, a measure of the temper- 
ature and pressure distributions can be obtained by observations around the 



GO till/, oxygen band. In each case, the quantum pinnies of the line formation 
is sufficiently well understood that the uncertainties in the atmospheric |«ra- 
meters dominate the uncertainties in the transition rates over most of the 
range of conditions of interest to us. 

In Figure I, we give the attenuation spectrum for the whole atmosphere to il- 
lustrate the contributions of the various com|x>nents of the atmosphere. Note 
that the residual IUO attenuation is due to the is.'t.ti ('ill/, water line. Further, 
over the region of interest to us <1.4 to 10 Cill/I the frequency dependence of 
the absorption is small (in a cloud free easel. In Figure 2. we give the bright- 
ness temperature spectrum corresponding to Figure 1. Note that the atmos- 
phere is opaque within the GO Gil/ oxygen band due to pressure broadening. 

The 22. 2.'tf> GHz water vapor line is still partiullv transparent, however. In 
order to observe the contributions from the higher altitudes, it is necessary 
to make use of the line wings particularly for the oxygen band. 

The time dependence of the atmospheric index of refraction, and therefore the 
path length correction, is dominated by the variations in atmospheric water 
vapor content. The variations in the* atmospheric water vapor content can Ik* 
thought of as corresponding to three levels of spaeial structure: 

1. Macrostructure: Tvpical sizes are of the order kilometers. These 
structures correspond to meteorological phenomena such as "ber- 
muda" highs and large cyclones. 

2. Mezzostructure: The characteristic size of these structures is a few 
hundred meters. These structures correspond to typical clouds and 
are the principal cause of atmosphere -induced fringe phase fluctua- 
tions with short -baseline (1-5 km) interferometer svstems. 

3. Microstructure: These structures are a few meters or less across 
and arc the cause of non-ionosphcric scintillations. It should be 
pointed out that the division into size regimes is strictly arbitrary. 
Except for stable cells, a continuum of sizes is almost always 
present. 

Only a relatively limited body of information has appeared on the characteris- 
tics of the various size scales. This is particularly true for 'he microstruc- 
ture. It is important that additional observation-* be obtained so that the de- 
tailed importance of each effect be determined. 

Current remote sensing svstems directed toward satellite applications are 
concerned mostly with the macroscale water vapor structure. Various 
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Attenuation Spectrum for the Entire Atmosphere under Cloud-Free Conditions 
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experiments < i . «*. , Itoscnkranz ct al , 1972: Ntaelln ct al, 1973; Wilhcit c*t al , 
19751 show that the variation of prcctpitublc water between two macrostruc- 
turea, while typically 12 to l-l 2 |(nt /cm : , cun la* an great an 4 gm 'em ’ 

The time scale for thin change de| tends on the S|tccd of th«* struct*. c's move- 
ment and cun In* as abort as 20 minutes or as long as several hours. Path 
length changes of nearly 1 2 meter ean occur in ihis period. The fronts which 
bound th« high and low pressurt regions often h ivc a layered structure due to 
the intrusion of warm air into cold air and \ Ice versa. As a result, the region 
oi the boundary between macrostructures in often a region of rapid and loca 
lion dependent path length changes. 

Short baseline interferometer observations provide the liest data on the mezzo- 
structure. it is possible, lor example, to determine both the size distribution 
and motion of water vapor ’'cells'* by observing the time correlation of phase 
fluctuations for several baselines oi a iven interferometer (Hinder, 1970; 
llasart rt al, 1170; Wessoling et al, 10i «». I*he»e studies show that the me*- 
zostructure causes path length variuti m >1 between 0. 1 cm and 5 cm and the 
characteristic size is In tween 300 and 1200 meters with a median size of 700 
meters (Hinder and It le, 197 ti. Since such structure is observed to occur at 
1 - 2 km and seems to me. e with the winds at altitude, we need to make meas- 
urements with a time seal*' as short as 5 seconds to distinguish the fastest 
moving of the cells. 

The microstructure is responsible for tropospheric scintillations. Virtually 
nothing has appeared regarding the nature of such structure but it is .ikclv 
that microstructure represents local Index el refraction variations within a 
cloud and perhaps corresponds to turbulent eddies. It would be difficult to 
separate the ionosphere and tropospheric scintillations in the wavelength range 
of interest to us since the two species of ii regularities w ill vield similar angu- 
lar fluctuation spectra. I'roquencv dependences would be required to affect a 
separation. 

The p roc dure foi determining the path length correctloi from (Missive micro- 
wave observations must of necessity Ik* a statistical one. This results from 
the lack of perfect coi relation between the distributions of meteorological vari- 
ables and tin measured brightness temperatures. dadiattve transfer calcula- 
tions show that individual frequencies have an altitude response that (teaks at 
an altitude which depends on the frequency distance from the two major line 
complexes. An example of these dependences, called weighting functions, is 
given in Figure Hen vvr plot the normalized weighting function: 

\V, (hi - or,, (hi ' 0 M>o <h> 
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Figure 3. Normalized Weighting Functions for Frequencies Around 

the 22.235 GHz Water Vapor Line 
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when* n is the Hluorplion confident as u function of altitude at frequency »• 

and p is the water vapor density as a function of altitude (after Staelin, 
If jt ) 

IDCfil. As can Ik* seen, the emission due* to atmospheric water vapor at any 
particular frequency conies mostly from a layer about 15 km thick. 


The inversion procedure must include all relevant processes and must treut 
the radlali u e Iransfei problem in detail. The principal confusing effect in the 
frequence rango of interest is the distribution ot liquid water. The liquid 
drops hi typical clouds mukc a very minor contribution to the path length since 
their si/e distribution contains no drops larger than about 50 mm. However, 
the contribution of the liquid water to the atmosphere brightness temperature 
ean In* over loo ot tin* liquid tree value. Therefore the path len;’*h correc- 
tion scheme must correct for the influence of the liquid water on the bright- 
ness temperature measurements. 


The inversion schenu is arranged so that the Influence ol radiometer noise 
and the lack ot pc* * correlation between the measured brightness tempera- 
tures and Hu* pa* rs of interest are both treated. The details of the sta- 

tistical InveiTwii procedure are given in Haul et ai (1972*. Briefly, the |iura- 
me*. rs of interest art obtained by a matrix operation: 


p T> • ti 


(5) 


whet.° p is the vecto. »l derived quantities, <1 ts a vector of linear combina- 
tions »l the observations and i> is the matrix which performs the inversion. 

The inversion matrix 1" in obtained by a statistical analysis ol the* computed 
brightness temperatures for n series of observed atmospheres whose proper- 
ties spun the range ol atmospheric conditions ot interest. The advantage of 
this approach is that the parameter vector p must lx* an Improvement over a 
parameter vector d* rived wholh from local statistics. The result of having 
noisy brightness temperatures or of having no measurements is to force p to- 
ward the local statistical value. Modifications of this simple scheme are re- 
quired to account for non-linearities and to avoid a singular f) but these modi- 
fications are based on sound statistical principles. 

An analysis ot this tyjK* can also be us *d as the basis of a simulation study to 
find the best instrument configuration and predict the ultimate accuracy ob- 
tainable. Schnper et at (1970) have performed such a study for a specific fre- 
quency complement (22, 24, >2. (»:j and ->:t. »> > tili7' and tound that, together 
with surface meteorological variables (humidity, temperature and pressure), 
a radiometer system with o. 2°K accuracy permits a determination of the zenith 


s 



path length correction to i. 1 cm. A similar study will be reported in Section 

ill. 


M. THE CURRELT STATE OF INSTRUMENT DEVELOPMENT FOR PASSIVE 
MICROWAVE RADIOMETRY 

Instruments employed in pussive rndiometry will seem insensitive bv the 
standards of present day radio astronomy. In *» radiometer, sensitivity is 
usually less important than absolute stability and accuracy; a useful radiometer 
for a remote sensing system must have exceedingly great stability and an ac- 
curacy (repr><ducibility) substantially better than most conventional radio as- 
tronomy receivers. Thus, an entirely different set of design criteria is re- 
quired for passive microwave radiometry. 

Front end amplifiers, such as parametric amplifiers, cannot be applied to 
radiometer receivers because their gains cannot be determined with sufficient 
accuracy and their absolute stability cannot be controlled. The traditional 
noire »ube calibration, even when attenuated to the same level as the incoming 
signal, is useless because the absolute ‘hermodynamic scale is not repeatable 
even when the noise tube remains on ccuinuo isly. Further, it is extremely 
difficult to determine the attenuation of a precision attenuator to the accuracy 
required for use with a noise tube calibration. 

For these anu other reasons, remote sensing radiometers are normally mixer 
receivers with Pi rt ke switching. latching ferrite circulators are used for 
switching and the entire radiometer is normally contained in a temperature 
controlled enclosure. Both the IMcke load and the enclosure are usually main- 
tained at a temperature slightly above ambient to improve thermal stability 
even though this may increase the power requirements. 

Instrumental calibration usually proceeds in two wavs. In the first procedure, 
the radiometer input (usually a horn antenna at the frequencies of interest tc 
us) is surrounded by an enclosure of Eccoscrb or similar microwave black 
body material of known temperature and the radiometer output is observed. 

This yields a calibration point at high brightness temperature (perhaps 280° k). 
The input is then immersed in a similar enclosure containing boiling liquid 
nitrogen. This yields a calibration point at low temperature (77°k). Tempeia- 
ture coefficients for the various waveguide and switch losses are obtained by 
observing the high brightness temperature load while the physical tempera- 
tures of the components vary. These measurements serve to set the absolute 
levels of the components independent of measurements of the sky brightness 
temperature. 
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To handle variations in the effective gain during observations, the* tempera- 
tures of all the temperature dependent components in front of the Dicke switch 
are monitored by thermistors. Additionally, a latching circulator Is located 
immediately alt of the antenna so that calibration loads can lie switched into 
the signal path in pluee <>1 the antenna. In modern receivers, these loads are 
always thermally compensated resistors. In the many systems, one is at a tem- 
perature higher than the bight nt brightness temperature expected while the 
other is lower than the lowest expected bright less temperature. Should it 
prove impossible to obtain a low calibration low enough, it is necessary to de- 
pend on a linear receiver characteristic to extrapolate two high calibrations 
downwa rd. 

In the past, satellite instruments have used either a radiation cooled Mad or 
else a horn viewing "cold" spuct to obtain a low reference. Ground based and 
aircraft instruments have used a boiling liquid nitrogen reference. A precision 
load and platinum resistance thermometer are immersed in a Dewar containing 
the boiling liquid nitrogen. Boiling liquid nitrogen has the advantage that the 
temperature is set by a phase transition whose temperature is very stable as 
Ion , as atmospheric conditions (pressure) are reasonable. Thus thermal gra- 
dients are suppressed, further, the thermal capacity of liquid nitrogen is not 
so high as to be a hazard unlens one insists i a being aggressively stupid. " 

Both varieties of •eferei loads must be careful I \ temperature controlled so 
that their emissivitv is in a riant and their temperatures are well defined. In 
practice, the best systems allow gradients in the loads of 0. 1 k at most. Th*J 
has proved to require considerable care but has yielded the advantage that the 
cold load Dewar need not be filled for nearly 5 hours. 

In the rest of this section, we will discuss the characteristics of two instru- 
ments which illustrate these design considerations. There are the microwave 
spectrometer now in operation on Nimbus 5 (NEMS) and the 37 GHz skv tem - 
perature radiomett r used on the NASA CV 990 during the Bering Sea 
Experiment. 

The Nimbus 5 microwave spectrometer (NEMS) consists of five 250 MHz band- 
width receivers centered at 22.235 GHz, 31.4 GHz, "3.05 GHz, 54.90 GHz and 
58.30 GHz. The purpose of the experiment is to sound the temperature, pres- 
sure and w..ter content of the atmosphere in nearly all weather conditions 
(Staelin et al , 1972; Staelin et al . 1973). Each receiver is a load -switched, 
double -sideband, superheterodyne receiver with a relative accuracy of 0. 1 - 
0.2^. In Figure 4, we give a block diagram of a tvpical NEMS receiver. The 
two calibration loads ere maintained at 215 k and 290 k while the Dicke refer- 
ence load is maintained a' ambient temperature. The instrumental time 
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Figure 4. Block Diagram of a typical NEMS Receiver 





constant is about 2 seconds and the calibrations are observed every 2. r »<> sec- 
onds. The horn antennas used at each frequency yield a beam width of 10°. 

The CV 5)90 skv temperature radiometer was used to determine the liquid 
water content of clouds (Wilheit et al . 197."») and for correcting surface emis- 
sion measurements for reflected sky emission (Webster et uj, 1975). A block 
diagram is given in Figure 5. The cold load is in boiling liquid nitrogen, 
while the hot load is maintained at about 30t)’k. The rms uncertainty in the 
measured skv temiiersture is 0.9°k for a 1 second integration. Because of the 
need for the highest possible absolute accuracy, the instrument only views the 
sky for 1/3 of the time. The remaining time is spent viewing the hot and cold 
loads. It should l>e pointed out that the cold load has performed without any 
difficulty even during short high-g turns and zero-g maneuvers. 

111. A SIMULATION STUDY OF THE PPME REQUIREMENTS FOR PATH 
LENGTH M EASU H EM ENTS 

Using the techniques outlined in Section I, we have simulated the |K*rformance 
of a path length measurement system and so optimized the choice of frequen- 
cies for the final sensor package. This study is based on a formulation devel- 
oped by T. T. Wilheit and T. C. Chang for the Seasat Phase A Study (1973) and 
the reader is referred to the Phase A report for details. 

To summarize, the atmosphere is divided into 100 lavers. The sensor |K»ck- 
age is assumed to include surface temperature, pressure and humicity meas- 
urements as well as microwave radiometers. The microwave sky tempera- 
tures were calculated from the best current transition probabilities and line 
formation theory. The atmospl.sre models were taken from the Handbook of 
Geophysics and Space Enviionments and span the range from arctic to tropical 
conditions. Cloud models were selected to span the range of most cloud types. 

Before reviewing the results, it is appropriate to summarize the implicit and 
explicit assumptions: 

1. The Handbook atmospheres are assumed to span the conditions of in- 
terest. Thus temperature inversions and other deviations from the 
models art not treated. In addition, temporal variations were not 
treated. 

2. Cloud attenuations have been assumed to depend only on the total liq- 
uid water . orient and Rayleigh scattering has been assumed. 
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Figure 5. Block Diagram of the CV 990 Skv Temperature Radiometer 












3. The wot ami dry components have not Iwen separated. The integral 
has been approximated by a sum. 

4. The radiometers have lieen assigned an rms measurement noise of 
0. 5 k. The accuracy requirements for the meteorological sensors 
will be discussed below. 

5. It has been assumed that errors in the line formation theory make no 
contribution 

We first examine the results for zenith observations. The first major conclu- 
sion is that the most obvious frequency set 422. GHz, either 18 or 31.4 GHz 
and 53.29 Gtlzi together with surface meteorology- measurements (humidity, 
pressure, temperature) determine the correction with an rms accuracy of 1.8 
cm. In agreement with Wilheit and Chang (1973) and Schaper et al (1970), we 
find that, it the 53. 29 GHz channel is not included, the error increases to 2.2 
cm. 


As t. test of the necessity of lower frequency channels, the H.4 GHz off-source 
measurements from the pronosod PPME continuum receiver were included 
with the I'-, 22.235, 53.29 triplet. If the sky brightness temperature is meas- 
ured with an rms noise of .). 1 k, the accuracy improves by a small amount 
(0. 1 ctn). Alternatively, if the more conventional rms noise of 0. 5°k occurs, 
the error increases by 0.2 cm. Therefore, unless the skv brightness temper- 
a‘ure is measured to the accuracy level of the best current radiometers (i.e. , 
the NKMS experiment; t ee Section II), the 8.4 GHz channel does not signifi- 
cantly improve the determination. 

An iterative scheme was used to select the best pair of channels near the 
22. 235 GHz water vapor line. The scheme shows that, regardless of the sec- 
ond frequency, the water vapor line frequency should be about 1 GHz from the 
line center, i e. . eith *r 21 or 23 GHz. yin improvement of o. 2 - 0.3 cm re- 
sults dope. ding on the second frequency. This Is as expo ted since the per- 
centage contribution to the total brightness temperature by higher altitude wa- 
ter vapor is greater ofl the iine center. The region about 1 GHz from the line 
center is an effective c ompromise between the need to measure the high alti- 
tude water vapor and the need to have a high signal -to-noise ratio. 

The second frequency is required to be at least 2 GHz away from the* near-line 
frequency. Since the second frequency gives the correction for the presence of 
liquid water, within limits, the farther it is from the near-line channel, the 
lower the error is. The 1 in' it is set on the low frequency side by the need for 
a high signal -t^ noise rctio and is about 17 GHz. I he limit is set on the high 
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frequency Hide by the increasing contribution from the oxygen hand at 60 (illt, 
and is about 40 GHz. 

In Table 1, we give the rms error in the path length determinations and the 
total correlation between the measurements and the path length for a selection 
of the frequency pairs which show the lowest errors. Note that the IS and 19 
GHz pairs * re essentially identical. An extension to 31 GHz results in a slight 
improvement. Because of the sharpness of the local minimum in the sky 
brightness temperature between the water vapor line and the oxygen, moving 
1 GHz or more away from :tl GHz notably degrades the performance. 


Table 1 

Accuracy of Path length Determination and Correlation of Measurements 


F requency 

Pair 

RMS Error 

Correlation 

(GHz) 

(Cm) 

(%> 

18 

21 

2.23 

97.5 


19 

21 

2. 19 

97.fi 


25 

21 

2.20 

97. 5 


18 

23 

2.20 

97.6 


19 

23 

2. 19 

97.6 


25 

23 

2.20 

97.6 


31 

21 

2. 19 

97.6 


31 

23 

2. 19 

97.6 



The above calculations have assumed that along with the microwave measure- 
ments, the surface humidity is measured to 5%, the rurface pressure to 1 mb 
and the surface air temperature to 2°K. Sincf the cos. of surface meteorolog- 
ical instruments is a very strong function of the desired measurement accu- 
racy, the contribution of improved accuracy for the meteorological instruments 
must be determined. To illustrate the change in performance for improved ac- 
curacy in the instruments, we give in Table 2 accuracies and correlations for 
the 18 GHz - 21 GHz pair as the accuracy of the meteorology instruments im- 
proves. Table 2 shows that the ultimate accuracy contribution of the meteor- 
ology instruments increases much more slowly with instrumental accuracy 
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Table 2 


lnl lu« nee *>l Surface Meteorology Instruments on Accuracy 


— — ■ ■ - 1 

Assumed Meteorology 
Error 

Path Length 
Error 

Correlation 

Per;ent 

Improvement 

llunuiliiv 

.V 

2.23 

cm 

97.4% 


Air Tempera* are 

i K 





Pressure 

2n - .b 





Hum iditv 

31 

2.04 

cm 

97.4% 

81 

Air Temperature 

! a K 





Pressure 

lmb 





Humidity 

n 

2. 0 cm 

97.4% 

! ioi 

Air Temperature 

1 ■’K 





Pressure 

0. 5mb 





Humidity 

2<- 

1.9H 

cm 

97.4% 

in 

Air Temperature 

. 0 1 K 





Pressure 

.Oil’ll) 





1 — 


than the cost. For example, according to T. A. Clark (personal communica- 
tion!, the cost ot the 'ust set of Instruments is a factor of nearly ten higher than 
the cost of th«“ third ret. Accordingly it is not cost-effective to obtain the ac- 
curacy level of the last set in the table. It is however desirable to obtain the 
accuracy level of the third set since this yields an improvement of 10'? for an 
estimated cost increme it of about a factor of 2 (T. A. Clark, personal 
communication) 

We next address the o.f- zenith problem. The calculations show that the path 
length varies with the secant of the zenith angle to high accuracy. This is as 
expected since the puih length is the equivalent of the optical air mass which 
also varies as the secant of the zenith angle. • ic inversion study shows that 
the error also scales roughly as the secant, lr Table 3 we give the errors and 
correlations for some ot the frequency pairs of Table 1 for zenith angles of 45 
and 10 . 

Clearly, by 60 zenith angle, the error has reached an unacceptably high level. 
The PPMK goal calls foi an accuracy of no worse than 4 cm total at large 
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Table 3 


Behavior of the Frequency Pairs for Off-Zenith Observations 


Frequency Pair 

/ 1 5 

Z « 60 ° 

(GHz) 

Error 

(Cm) 

Correlation 

<%> 

Error 

(Cn»> 

Correlation 

<%) 

18-21 

4.22 

95.2 

6.00 

95 . 1 

19-21 

4.24 

95 . 1 

6. 11 

94.9 

25-21 

4.39 

94.8 

6.21 

94.8 

18-23 

4.28 

95.0 

6.24 

94.7 

19 —?' i 

4.34 

94.9 

6.43 

94.4 

25-23 

4.46 

94.6 

6 . 50 

94.2 

31-21 

4.21 

95 . 1 

5. 99 

95 . 1 

31-23 

4. 14 

95.3 

6.28 

94.6 


zenith angles. The question arises as to how much of th" error is due to the 
wide span of atmospheric conditions included in the analysis, in the preceding 
analysis, the total range of atmosphere models covers conditions from clear 
polar winter to humid cloudy tropic summer. We now restrict the range of at- 
mosphere models to mid-latitude conditions. For each case in Table 1 and 3, 
the errors are cut by a factor of 2. To illustrate the improvement, we have 
calculated the 19 and 21 GHz brightness temperatures and path length for the 
U.S. standard atmosphere and performed the inversion. Since the U.S. stand- 
ard atmosphere represents the average U.S. conditions, we ex|>ect that the in- 
version will recover the calculated path length to high accuracy. As Table 4 
shows, this is the case even for heavy clouds. Therefore, it may be necessary 
to develop the regression constants for each of the PPME sites fron leasured 
atmospheric conditions which span the total range of conditions at the site. 

Analyses were also made for three frequency schemes to indicate the influence 
of the increased path length as the zenith angle rises. The two lower frequen- 
cies were fixed at 18 GHz and 21 GHz while the third frequency ranged from 37 
to GO GHz. The region from about 55 GHz to 60 GHz made no contribution to 
the Inversion since the atmosphere is opaque within the oxygen band. In the 
region from 37 to 50 GHz, the contribution is of the same order as at 31 GHz 
and becomes progressively worse as the frequency rises. In the region between 
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Tabli- 4 


Inversion til U.S. Standard Atmosphere 


19 (ills, 21 CHz, Surface Pressure, Iluinldttv, Air Temperature 


— — r 

z = 0 

0. 1 ent No Clouds 

j 

0. 1 cm Heavy Clouds 

Z = 45° 1 

«). G cm No Clouds 

0. 7 cm Heavy Clouds 


50 and 55 GHz, the brightness temperature is proportional to the temperature 
and pressure distributions. Therefore, the error decreases somewliat for a 
thiid frequency in this range. As the zenith angle increases, the correlation 
of the third frequency decreases. This is due to the atmosphere becoming op- 
tically thick at the larger zenith angles between •»<-» and 55 GHz. To illustrate 
this, we give in Table 5 the errors including correlations for the l* Gils, 21 
GHz, and 59 GHz triplet for zenith, 45 and GO zenith angle. 


Table 5 

Error and Individual Correlations for the IN GHz, 21 GHz and 53 GHz Triplet 


Zenith Angle 

Error 


Correlations (%) 


(Degrees) 

(Cm) 

IH GHz 


21 GHz 

<■!!/ 

0° 

1.8 Cm 

84% 


cn 

83% 

45° 

3.7 

84% 


57% 

50% 

60° 

5.8 

84% 


57% 

20% 


Thus, a system which will measure the path length correction to the required 
accuracy should be configured as follows: 

1. At least two microwave radiometers (IN GHz, 21 GHz) with 0.5 K rms 
noise and 1 J K long term stability. If a third frequency can be accom- 
modated it should be about 53 GHz. 

2. Surface meteorology instruments which measure relative humidity to 
3 , air temperature 1 K and pressure to 1.0 mb. 


IN 
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